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H '■ 

C^ , We use the galaxy angular power spectrum at 2 ~ 0.5 — 1.2 from the Canada-France-Hawaii- 

Telescope Legacy Survey Wide fields (CFHTLS- Wide) to constrain separately the total neutrino 
mass Yl i^<^ fi'iid the effective number of neutrino species A^off ■ This survey has recently benefited from 
^S| I an accurate calibration of the redshift distribution, allowing new measurements of the (non-linear) 

^Nj , matter power spectrum in a unique range of scales and redshifts sensitive to neutrino free streaming. 

Our analysis makes use of a recent model for the effect of neutrinos on the weakly non-linear matter 
^^ power spectrum derived from accurate N-body simulations. We show that CFHTLS, combined 

Q_^ , with WMAP7 and a prior on the Hubble constant provides an upper limit of X^m^ < 0.29 eV 

r) ■ and A^eff = 4.17tj'26 (2^'" confidence levels). If we omit smaller scales which may be affected by 

non-linearities, these constraints become ^ nii, < 0.41 eV and Ncs = 3.981^'2q (2 a confidence 
levels). Finally we show that the addition of other large scale structures probes can further improve 
these constraints, demonstrating that high redshift large volumes surveys such as CFHTLS are 
O ' complementary to other cosmological probes of the neutrino mass. 

C3 ■ I. INTRODUCTION 
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Determining the neutrino mass is one of the great unsolved problems of modern particle physics. The standard 
^ ■ model contains three massless neutrino species; observations of neutrino oscillations in atmospheric and solar neutrino 
experiments have confirmed that neutrinos are massive, but cannot pin down their absolute mass scale. Fortunately, 
the Universe offers a new laboratory for neutrino physics. Massive neutrinos affect both the background expansion 
(/— , I and the growth of cosmological structure, making cosmological observations an unrivaled probe of the total neutrino 

• , mass fll 14I . 
L«; ' The redshift at which massive neutrinos become non-relativistic alters the time since matter-radiation equality, thus 

-vj ,. changing the position of the Cosmic Microwave Background (CMB) anisotropy peaks. Measurements of the CMB 
y—{ • have used this effect to find X^'^i' < 1-3 eV (95% C.L.) [j-Q' ^^^ ^ sensitivity of a{J2i^v) ~ 0.3 eV could soon be 
!• [ achieved with the Planck satellite (e.g. [4]). 
. ^ • Massive neutrinos also play a role in the formation of large scale structure. Once non-relativistic, they damp 

S^ ' the growth of perturbations within their free streaming scale, resulting in a suppression of the small-scale linear 
^ . matter power spectrum of AP/P ^ —Sfl^/flrn H- Measurements of the matter power spectrum can thus improve 
constraints on the neutrino mass considerably. Many analyses have been performed combining CMB data with Large 
Scale Structure (LSS) probes such as the 2dF Galaxy Redshift Survey (2dFGRS) [l-O, the Sloan Digital Sky Survey 
(SDSS) [ill-Ill, WiggleZ [H and the SDSS Lyman-a forest [H-il to constrain Y. m„. The combination of the SDSS 
DR8 LRG angular power spectra, WMAP7 data and an HST prior on the Hubble constant gives ^ nii, < 0.26 eV 
(95% C.L.), assuming a fiat ACDM model with the standard model effective number of neutrino species, Nes = 3.04 

A detection of iVcff > 3.04, as recently hinted at by the Atacama Cosmology Telescope (ACT) and the South Pole 
Telescope (SPT) experiments |3ll. |32||. would imply additional relativistic relics or non-standard neutrino properties. 
The additional energy density of extra relativistic particle species would change the redshift of matter-radiation 
equality, leaving imprints on the CMB anisotropics and matter power spectrum. Parameter degeneracies limit the 
constraining power of the primary CMB anisotropics on Ngg 0, [^, |3J|. The 95% lower limit is A^eff > 2.7 from 
WMAP7 alone 'u\. Adding information on the matter power spectrum and the Hubble constant can tighten the 
constraint significantly [2ll [23l [32. 34-37] . as can measurements of smaller scale CMB fiuctuations. 

In this work, we introduce new constraints from measurements of the angular power spectrum of galaxies at 



z ~ 0.5 - 1.2 in the Canada-France-Hawaii- Telescope Legacy Survey Wide fields (CFHTLS-Wide) 'S^]. This is 
the deepest wide-field survey of its kind, covering 133 sqr. deg. and sampling a comoving volume of 0.2 Gpc'^ ' 
from z = 0.5 — 1.2. This data has been used for cosmological studies, including weak lensing constraints 
1421 . In particular, a joint analysis using the weak lensing measurements has given a limit for massive neutrinos of 
^m,y < 0.54 eV (95% C.L.) [4j|. The small-scale clustering was studied using the halo model by ref. [4J] and the 
deprojected power spectrum was presented in ref. ;455. 

Considerable gains can be made by probing the power spectrum on weakly non-linear scales, provided one can model 
the galaxy bias sufhciently accurately. On large scales, the galaxy power spectrum is well-fit by the dark matter power 
spectrum and a linear galaxy bias: 6g = bgSdm- On small scales, where the details of galaxy formation come into play, 



this relation breaks down and a general, scale-dependent, galaxy bias model is necessary [46l - l48l | . This results in a 
degeneracy between galaxy bias and cosmological parameters that is difficult to break without further observations 
such as gravitational lensing |49| . 

The scale-dependence of galaxy clustering has been investigated with simulations [4^ [501 ^^^ observations ^31 ■ 
These studies show that galaxies with lower luminosity have a weaker scale-dependence of the bias than the most 
luminous ones at fc > O.lh/Mpc. Additionally, the scale dependence becomes more severe for strongly biased tracers; 
negligible scale dependence of the halo bias is found in simulations at k = 0.15 h/Mpc for halos with bg < 2, while for 
more massive halos with bg > 2, there is a 20% effect [50|. For the SDSS main sample, the Halof it non-linear matter 
power spectrum with a constant bias factor has been demonstrated to be a good fit to the galaxy power spectrum 
|5l| . while the power spectrum of the LRG sample diverges at k = 0.2h/Mpc |48| . 

We select the range of scales to study based on the dominance of the two-halo term in the halo occupation dis- 
tribution (HOD) model of the power spectrum. In ref. [4J], HOD fits are made to the correlation function of the 
CFHTLS galaxies. We have checked explicitly that the two-halo term is greater than the one-halo term to £ = 960 
in our lowest redshift sample (0.5 < Zphot < 0.6). In Sec. |V]we test our CFHTLS data set using a comparison of the 
constant bias model with a two-parameter model in the context of ACDM without massive neutrinos. We conclude 
that a constant bias model is sufficient. 

The CFHTLS data set is ideal for this study because (i) at higher redshift the onset of non-linear growth happens 
at smaller scales; and (ii) we are targeting "normal" galaxies: the flux limited sample that we construct (zab < 22.5) 
selects galaxies to z = 1 with luminosities of Mg ^ —20, similar to the SDSS main sample at z = 0.1 [52ll53j. 

Since we analyze the projected density field, any modulation of the power spectrum due to redshift-space distortions 
is minimized, although these effects must be considered for surveys covering a larger fraction of the sky than CFHTLS 
[19| . We must also consider systematic errors arising from the luminosity dependence of the galaxy bias. In a flux- 
limited sample with a mixture of galaxy types, luminosity-dependent biasing can modify the slope of the galaxy power 
spectrum. However, for current surveys such as 2dFGRS, the effects can be neglected |54| . 

In this paper we will present constraints on the total mass of neutrinos, ^ m^, and the effective number of neutrino 
species, N^s, from different cosmological observations. The structure of the paper is as follows: In section HIl we 
briefly review the general formalism of the angular power spectrum C^. In sections Hill and HVl we describe the recently 
released CFHTLS clustering power spectra data set, and the other data sets we use. Section |V] contains our main 
results on J^''^'^ ^^-^ Ncs, while section IVTl is dedicated to the conclusions and discussion. 



II. ANGULAR POWER SPECTRUM 

The angular power spectrum, Ci, is a projection of the spatial power spectrum of fluctuations, -Pgai(fc, z), where k 
is the comoving wave number and z is the redshift. The equation for the projection is 



26^ r 
Ct = / k^dkP(k)g\(k) 



(1) 



where P(A:) is the matter power spectrum today and 6 is the assumed scale-independent bias factor relating the galaxy 
overdensity to the mass overdensity. In our calculations, we assume the values of bias are constant in each redshift 
bin; with width Az ~ 0.2, these are relatively narrow. We also neglect the effect of redshift space distortions, as their 
effect is limited to large scales, £ < 30 [53|, and our largest bin is at ^ = 40. With these assumptions, the kernel giik) 
is given by 
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where je{x) is the spherical Bessel function, x{z) is the comoving distance to redshift z, and dN/dz{z) is the normalized 
redshift distribution of the survey. 



The Limber approximation is valid for all but the largest angular scales {i > 10). For £ ;:^ 1, we have 

^ J edkMkx)ji{kx') = ^Six ~ x') , (3) 



and the equation for the projection becomes 

where dVc{z) is the comoving volume element, dVc{z) — x^dx/dz 



dz 



'^' (4) 
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III. Ci MEASUREMENTS OF CFHTLS 

A. Data 

We analyze the galaxy density maps constructed by ref. ^] from the Canada-France-Hawaii Legacy Survey Wide 
fields (CFHTLS- Wide). The data set is based on the CFHTLS T0006 release, including photometric redshift estimates 
[531 • A selection of galaxies was made photometrically to an apparent flux limit of iab = 22.5 in three photometric 
redshift bins: 0.5 — 0.6, 0.6 — 0.8, and 0.8 — 1.0, labeled S6, S7 and S8. We are not overly concerned with the error in 



the photometric redshift estimates, although it is typically ^z/{l -\- z) = 0.04 [52|; we only require knowledge of the 
redshift distributions of the samples. In ref. [45|, the redshift distributions were measured using a subset of ^^14000 
spectroscopic redshifts obtained from the VIMOS Public Extragalactic Redshift Survey (VIPERS) [5g|. We also use 
those redshift distributions here. 

B. Cii estimator 

We estimate the C^ with a maximum likelihood approach first applied to cosmic microwave background maps 
[571 l5S| . It is well suited to fields that may be described by a Gaussian likelihood function, and it has been extended 
successfully to measure the angular power spectrum in galaxy surveys [59, ^] . 

For surveys that cover only a relatively small patch of sky, such as CFHTLS, the maximum likelihood approach has 
advantages over other techniques. Generally, it is challenging to measure harmonic modes that are strongly affected 
by the survey geometry. Modes become correlated due to the convolution effect of the survey mask; further, the 
hard edges can introduce 'ringing' in harmonic space. Our estimator accounts for these effects and measures the Ci 
with minimum variance (in the case of a Gaussian field). Along with the C^, we derive the covariance matrix of the 
measurements and window functions. 

The galaxy density maps were constructed on a HEALPix^ grid with an angular resolution of 7 arcmin (ugido = 512) 
|6l| . We order the m pixels of the density map, form a data vector, x = [5(no),5(ni), ..., (5(n,„_i)], and write the 
covariance of the data as Cy — (xiXj). The pixel covariance matrix is given by the sum of the signal and the noise 
components 

97 I 1 

Cy = ^ ^— T'Kcos 9,j)BfCi + N,, , (5) 

where Vi are the Legendre polynomials, Oij is the separation between pixels i and j and Nij is the noise covariance 
matrix. The noise matrix is taken to be diagonal with Poisson elements given by Nu = wf /n, where Wi is a weight 
accounting for partially-sampled pixels. The finite resolution of the pixelised map attenuates the power spectrum by 
the pixel window function, B^, which depends on the pixel geometry [6ll|. 

To simplify the following expressions, we introduce a symmetric matrix P;,. This is a sum of the Legendre polyno- 
mials over a band-power indexed by h including [ff,,£b+i)- The components are 

Pb,^J= E ^^^KcOS%)i?^ (6) 

lb<l<lb+i 



^ |http : //healpix . jpl . nasa ■ gov] 



TABLE I: CFHTLS samples 





Scale ( h/Mpc) 


Photo-z sample z 


^ = 630 £ = 960 


S6 


0.5 < zphot < 0.6 


0.557 


0.43 0.65 


S7 


0.6 < Zphot < 0.8 


0.687 


0.36 0.55 


S8 


0.8 < Zphot < 1.0 


0.839 


0.31 0.47 



The quadratic band-power estimator is given by 

Cb^\Yl ^-J {x^C-iPfeC-ix - Tr (C-iPfeC-iN) } 



(7) 



The matrix A is a mixing matrix that sets the normahzation and may be specified to form linear combinations of the 
bin estimates. It is related to the Fisher matrix, 



Fw = ixr (C-iPfcC-iPfc,) 



The expectation and variance of the estimator are 



(A) = AFA , 



(8) 



(9) 



Var(A, A) = AFA^ 



(10) 



and the window functions are W — AF. The matrix A can be chosen in a variety of ways to optimize the estimator 
[57| . Because the Fisher matrix is ill-conditioned due to the small survey size, we use the robust normalization 

The covariance between the C^ measurements of two redshift slices labeled A and B may be estimated as 



Var(Cf,Cf) 
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where Cf^ is given by eq. ([T]), rewriting as 



C/^ = ^^^ / k'dkPik)gfik)gfik) . 



(11) 



(12) 



We must further convolve by the survey window functions. This is an idealization that neglects the precise survey 
geometry but we may scale it to match the variance computed from the quadratic estimator 



Ct measurements 



The quadratic Ce estimator operates under the assumption that the galaxy density can be described by a Gaussian 
random field, which is not valid on small scales. However, Gaussianity can be a reasonable model for galaxy counts 
in projection even to non-linear scales. Our map resolution is 7 arcmin/pixel (ugide = 512) and, in the first redshift 
bin, S6, the mean and median counts are 33.6 and 32.0 galaxies/pixel, respectively, with skewness 1.18. At a higher 
resolution of 3.5 arcmin/pixel, the mean and median become 8.42 and 7.72 with skewness 1.48. We see that the 
map resolution of 7 arcmin/pixel gives a good balance between Gaussianity and angular scale. We expect that the 
Gaussianity assumption is still valid at scales typical of the pixel size at £ = 1500. 

We compute the angular power spectrum in 48 bins over the range i = 2 — 960, with width Ai = 20. The Fisher 
matrix and window functions were computed with step size Ai =10. Due to the small angular extent of the survey, 
we cannot probe low £ and, indeed, the computed window functions show that we are not sensitive to ^ < 40. 

We estimate Ci for the three photometric samples, S6, S7 and S8. The measurements are shown in figure [TJ The 
Fourier scales that are probed in the angular power spectrum depend on the projection kernel gg. In figure [21 we plot 



representative kernels. We use two limits in angular scale in the subsequent cosmological analysis: ^j, 

-^max — you. 
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FIG. 1: The angular power spectra for the three different redshift bins. We show the data with la error bars and the linear 
(dashed curves) and non-linear (continuous curves) theoretical angular power spectra. 



IV. EXTERNAL DATA SETS 

Besides the angular power spectra Ce of the CFHTLS measurement, we will also consider the following cosmological 
probes: i) power spectra of CMB temperature and polarization anisotropics; ii) power spectra of luminous red galaxies; 
iii) measurement of the current Hubble constant; iv) luminosity distances of type la supernovae. These data sets are 
all well established cosmological probes that have been extensively investigated and already provide tight constraints 
on the cosmological parameters of the concordance ACDM model we will use here. 



A. CMB Power Spectra Data 



To incorporate the seven-year WMAP (WMAP7) CMB temperature and polarization power spectra, we use the 
routines for computing the likelihood supplied by the WMAP team [6|. The WMAP7 polarization data are composed 
of TE/EE/BB power spectra on large scales {2 < £ < 23) and TE power spectra on small scales (24 < £ < 800), while 
the WMAP7 temperature data includes the CMB anisotropics on scales 2 < £ < 1200. 

Here we do not use other small-scale CMB temperature power spectra measurements, as adding them would not 
significantly improve the constraints on the cosmological parameters, especially those on the total neutrino mass 
focussed on in this paper. 
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FIG. 2: The projection kernels ge{k) defined in eq. ([2]) for tlie tfiree CFHTLS samples S6, S7 and S8 with mean redshifts 0.557, 
0.687 and 0.839. The normalized kernels are plotted. 



B. Power Spectrum of Luminous Red Galaixies 



The power spectrum of LRGs measured by SDSS is a powerful probe of the total mass of neutrinos, '^m^, and 
the effective number of neutrino species, A^cff- We thus include the SDSS DR4 LRG power spectrum [6^] which has 
a mean redshift z ^ 0.35 and use data points on scales 0.012 /i/Mpc < fccff < 0.203 fe/Mpc in the analysis. 

We checked the constraining power of the SDSS DR7 LRG power spectrum I2l| and found that the constraint 
on the total neutrino mass does not improve significantly. For simplicity, we still use the SDSS DR4 LRG power 
spectrum [G^. Furthermore, we do not include the BAO information i22], as the measurement of BAO and LRGs 
power spectrum cannot be treated as independent data sets. 



C. Hubble Constant 

In our analysis, we add a Gaussian prior on the current Hubble constant given by ref. [63]; Hq = 74.2 ± 
3.6kms~^ Mpc~^ (68% C.L.). The quoted error includes both statistical and systematic errors. This measure- 
ment of Hq is obtained from the magnitude-redshift relation of 240 low-z Type la supernovae at z < 0.1 by the Near 
Infrared Camera and Multi-Object Spectrometer (NICMOS) Camera 2 of the Hubble Space Telescope (HST). This is 
a significant improvement over the previous prior, Hq = 72 ± 8kms~^ Mpc~^, which is from the Hubble Key project 
final result [6J]. In addition, we impose a weak top-hat prior on the Hubble parameter: Ha e [40, 100] kms^^ Mpc^^. 



D. Luminosity Distances 



Finally, we include data from Type la supernovae, which consists of luminosity distance measurements as a function 
of redshift; Di^{z). In this paper we use the latest SNIa data sets from the Supernova Cosmology Project, "Union 
Compilation 2.1", which consists of 580 samples and spans the redshift range < z < 1.55 [651). This data set 
also provides the covariance matrix of data with and without systematic errors. In order to be conservative, we use 
the covariance matrix with systematic errors. When calculating the likelihood from SNIa, we marginalize over the 
absolute magnitude M, which is a nuisance parameter, as done in refs. [66J, |67|. 




FIG. 3: Marginalized one-dimensional and two-dimensional likelihood (1, 2 a contours) constraints on the total neutrino mass 
and the three CFHTLS bias parameters from WMAP7-hHST-hSDSS-hSN-hCFHTLS data combination for ^^ax = 630. We also 
show the constraints on the bias parameters, assuming massless neutrinos (red dashed lines). 



V. NUMERICAL RESULTS 



We model the nonlinear matter power spectrum P{k) using the Halof it formulae [681] as modified to account for 
massive neutrinos 69]. This modified version of Halof it was obtained from an extensive suite of N-body simulations 
which treat massive neutrinos as an independent set of particles |7C| , and was particularly focused on exploring small 
scales at redshifts z = — 2. Similar simulations have been used to estimate the redshift space distortions and bias 
between matter and haloes [7l|. 

We perform a global fitting of cosmological parameters using the CosmoMC package [T^l , a Markov Chain Monte Carlo 
(MCMC) code. We assume purely adiabatic initial conditions and a flat ACDM Universe, with no tensor contribution 
to primordial fluctuations. The following six cosmological parameters are allowed to vary with top-hat priors: the 
dark matter energy density parameter f2c^^ £ [0.01,0.99], the baryon energy density parameter fib/i^ £ [0.005,0.1], 
the primordial spectral index rig £ [0.5, 1.5], the primordial amplitude log[IO^°As] G [2.7,4.0], the ratio (multiplied by 
100) of the sound horizon at decoupling to the angular diameter distance to the last scattering surface Gg G [0.5, 10], 
and the optical depth to reionization r e [0.01,0.8]. The pivot scale is set at fcgo = 0.05 Mpc~^. Besides these six 
basic cosmological parameters, we pay particularly attention to the neutrino mass fraction at the present day 
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the effective number of neutrino species, iVeff , and the three CFHTLS galaxy bias parameters, &S6, 6s' 
in our analyses we do not vary ^ ra^, and N^f^ simultaneously, since they are no longer degenerate 
assume Nes = 3.04 to constrain the total mass of neutrinos and ^ rrii, = to constrain Nes- 



(13) 

6s8. Note that 
. Instead, we 



A. CFHTLS galaxy bias 



We first measure the galaxy bias in the context of the ACDM model, without introducing massive neutrinos. We 
find the 68% C.L. constraints on the bias parameters to be: 635 = 1.41 ± 0.05, 637 = 1-54 ± 0.05 and 638 = 1.73 ± 0.06, 
which are in goo d agreement with those derived by ref. [45| . The bias values are consistent, but higher than those 
found by ref. [44| for volume limited samples in CFHTLS covering the same redshift ranges. This discrepancy could 



TABLE II: The 95% confidence level upper limits on the total mass of neutrinos 5I^?n„ from different data combinations. 



95% C.L. Y.m, [eV] 


Without HST Prior 


With HST Prior 




^max = 630 fmax = 960 


^max = 630 ^n,ax = 960 


WMAP7 


1.17 


0.50 


WMAP7 + CFHTLS 


0.64 0.43 


0.41 0.29 


WMAP7 + SDSS + CFHTLS 


0.47 0.35 


0.35 0.28 


WMAP7 + SDSS + SN + CFHTLS 


- 


0.33 0.27 



arise from differences in tfie analysis metliod; in ref. 4J], the cosmology was fixed and the fit was carried out over 
halo model parameters using the correlation function on small scales (< 1°). Generally, the best-fitting bias values 
are correlated with the other parameters considered in the analysis. In particular, we will find that they shift when 
we introduce massive neutrinos. 

We also check whether a scale-dependent bias is supported by the data, again assuming ACDM without massive 
neutrinos. We compute the best-fitting parameters assuming a scale-dependent Q bias model 74] with the form 



We fix the parameter A — 1.7Mpc/h and find the best-fitting Q. 

We first perform the fit using the linear power spectrum, such that P{k) = 6'^(A;)Pii„(fc). The best-fitting Q value is 
Qiin = 10.3 ±1.5. If we instead use a non-linear Halof it power spectrum (without massive neutrinos), the best-fitting 
Q value is reduced to Qhaiofit = 2.5 ± 1.2 and scale dependence in the bias at the < 5% level at fc < 0.6h/Mpc, which 
is negligible. The robust fit using Halof it and a simple constant bias is expected, since the data is well- fit by the 
two- halo clustering term on these scales [44>|. This demonstrates that the assumption of constant galaxy bias is a 
good one, and our results for massive neutrinos are robust. 

Before presenting the constraints on the total neutrino mass, we examine the degeneracies between the total 
neutrino mass and the three CFHTLS bias parameters. In figure [3] we show the marginalized one-dimensional and 
two-dimensional likelihood constraints from the WMAP7-I-HST-I-SDSS-I-SN-HCFHTLS data combination. The black 
solid lines and the red dashed lines denote the constraints with massless and nrassive neutrinos, respectively. When 
including massive neutrinos, the constraints on the bias parameters are obviously weakened and the median values 
are shifted: &S6 = 1.49 ±0.07, 6s7 = 1.62 ±0.08 and 5s8 = 1.81 ±0.08 (68% C.L.). This results in a strong correlation 
between ^ niy and the bias parameters; clearly shown in the two-dimensional contours of the last row of figure 
|3l A larger total neutrino mass will suppress the amplitude of the matter power spectrum on small scales further. 
Countering this while still matching the CFHTLS data on small scales requires larger bias parameters. Thus, further 
improving the constraints on the total neutrino mass will require better determination of the bias parameters. 

B. Neutrino Mass ^niv 

In this subsection, we present the 95% confidence level upper limits on the total mass of neutrinos from different 
data combinations after marginalizing over the other parameters, as shown in table HIl 

We start by presenting the limits using scales up to ^max = 630. Due to the strong degeneracies present between 
cosmological parameters, primary CMB anisotropics alone can place only weak constraints on the total neutrino mass. 
In the flat ACDM framework, WMAP7 alone constrains ^yn^, < 1.17 eV at the 95% confidence level Q. Adding 
the low redshift CFHTLS measurements breaks these degeneracies and the constraint on the total neutrino mass 
significantly improves to 

J2mu< 0.64 eV (95% C.L.) (15) 

for the combined WMAP7 and CFHTLS datasets. In figure 21 we plot the one-dimensional marginalized distributions 
on some cosmological parameters from different data combinations. As can be seen, adding CFHTLS data improves 
the constraints on the present matter density rim , the amplitude of fluctuations as and the bubble parameter Hq ; their 
68% confidence levels are shrunk from 0.334 ± 0.053, 0.711 ± 0.062 and 65.1 ± 3.7 for WMAP7 alone to 0.320 ± 0.043, 
0.759±0.044 and 66.3±3.1 for WMAP7+CFHTLS data, respectively. In figure[5]we show the two-dimensional contours 
in the (ctsiX) "^i')! (^m,X]"^i^) ^^"^ [Hvi^'m^) planes from different data combinations. Using WMAP7+CFHTLS 
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FIG. 4: One-dimensional marginalized likelihood constraints on the total neutrino mass ^ tti,^, as well as other cosmological pa- 
rameters, Qbh^, ^ch^, ^m, erg and Ho from different data combinations: WMAP7 alone (black solid lines), WMAP7+CFHTLS 
(blue dashed lines) and WMAP7-HSDSS+CFHTLS (green dotted lines) for £niax = 630. 



data (blue dashed lines) reduces the correlations between ^ m^ and other parameters and gives tighter constraints 
on '^mu than WMAP7 alone (black solid lines). 

Using the SDSS DR4 LRG power spectrum data [62], we obtain a constraint on the total neutrino mass of ^ rn^ < 
0.62 eV at the 95% confidence level, comparable to the limit from WMAP7-I-CFHTLS. For comparison, we also use 
the SDSS DR7 LRG power spectra data ^21] to constrain ^ m^ and get a 95% C.L. upper limit of ^ m^ < 0.61 eV, 
which is almost identical with that from DR4 LRG data. Due to this negligible improvement, in our calculations 
we still use the SDSS DR4 data to avoid using the more complicated DR7 likelihood code. Figure |4] (green dashed 
lines) clearly shows that the constraints on J7,„, ag and Hq are much tighter when including the SDSS data set. We 
obtain a constraint on the total neutrino mass of 'Y^rrii, < 0.47 eV at the 95% confidence level from the combination 
of WMAP7+SDSS-hCFHTLS data. 

We next consider the constraints using WMAP7 with an HST prior. The degeneracy between Hq and Yl "^t^; clearly 
shown in figure \E\ means that independent measurements of Hq produce significantly tighter overall constraints. 
WMAP7-HHST gives the 95% C.L. upper limit Y.'^'^ < 0.50 eV. After adding the low redshift CFHTLS clustering 
power spectra data, the constraint becomes even tighter 



^m^<0.41eV (95% C.L. 



(16) 



In figure [HI we show the one-dimensional marginalized distributions on some cosmological parameters from different 
data combinations when the HST prior is used. Clearly, the CFHTLS data improves the constraints on some cosmo- 
logical parameters, reducing the correlations between them and ^ m^ further, and thus significantly improving the 
limit on ^ m^. 

For comparison, we use the SDSS DR4 LRG power spectra data to constrain ^ rrii, and obtain the 95% upper limit 
0.45 eV, comparable to that from the CFHTLS data set. Similarly, we use the SDSS DR7 LRG data and obtain the 
constraint X) '^'^ < 0.43 eV (95% C.L.), which is consistent with previous work [6]. When using SDSS and CFHTLS 
data sets together, we get tighter constraints on the total neutrino mass; Y'^^i' < 0.35 eV at the 95% confidence 
level. Finally, we add the "Union 2.1 Compilation" supernovae data into the WMAP7+HST+SDSS+CFHTLS data 
combination and the constraint on the total neutrino mass becomes 



^m^<0.33eV (95% C.L.) 



(17) 



We now present the constraints obtained from extending the ranges of multipoles to smaller scales; up to t-n 
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FIG. 5: Two-dimensional contours in the (as,'^m,^), {Qrn,'^rn^) and {Ho,^m,^) panels from different data combinations: 
WMAP7 alone (black solid lines), WMAP7+CFHTLS (blue dashed lines) and WMAP7+SDSS+CFHTLS (green dotted lines) 

for ^max = 630. 



The improvements from these scales are not negligible: if we consider CFHTLS data in combination with WMAP7 we 
go from J2 "^i^ = 0-64 eV to ^ m^ < 0.43 eV. If we also add the HST data we obtain ^ niij < 0.29 eV, which has to 
be compared with the ^max = 630 results of ^ rrii, = 0.41 eV. Note that by adding SDSS and SN data the constraint 
improves and becomes '^ttIi, < 0.27 eV at the 2ct C.L. level, which is the tightest bound on neutrinos presented in 
this paper. 

It is evident that exploring the smaller scale matter power spectrum can significantly improve constraints obtained 
from larger scales. Our result above can be compared to the constraint obtained using only multipoles up to the linear 
theory value ^max = 150, which is ^ nii, < 0.43 eV for the combined data set. 

Finally, we have explored the effect that the error on the non-linear fitting formula for the matter power has on the 
final results. This error has been quantified in ref. 69], which presents an analytic expression, E{k,z), the error on 
the suppression in the power spectrum due to neutrinos, which depends on scale, redshift and neutrino mass fraction, 
and is at the level of 5%. We computed the constraints by considering P{k) — P{k)[l + aE{k, z)], with a being a free 
parameter with standard deviation of unity. For all data sets the 2 a upper limit for £niax = 960 weakens by only O.OleV, 
demonstrating that the numerical errors for the scales considered here are smaller than statistical uncertainties on 
the data point and will not impact the final constraints. Halof it will also underestimate the matter power spectrum, 
in a manner independent of neutrino mass, by around 5% [75|. Neglecting this error, because it underestimates the 
matter power spectrum, is conservative and will only produce slightly weaker constraints. Furthermore, since the 
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FIG. 6: One-dimensional marginalized likelihood constraints on the total neutrino mass '^mv, as well as other cosmological 
parameters, Q.^ and a% from different data combinations: WMAP7+HST (black solid lines), WMAP7+HST+CFHTLS (blue 
dashed lines), WMAP7+HST+SDSS+CFHTLS (green dotted lines) and WMAP7+HST+SDSS+SN+CFHTLS (red dash-dot 
lines) for fmax ~ 630. 



magnitude of this error is similar to i?(fc, z) and independent of neutrino mass, its effect should be small. 

TABLE III: The 68% and 95% confidence levels on the effective number of neutrino species, A'ctt from different data combina- 
tions. 



Data sets 


iVcff 




^max = 630 ^n,ax = 960 


WMAP7 alone 
WMAP7 + SDSS 
WMAP7 + CFHTLS 


>2.84 (95%C.L.) 

A 1o + 1.20/+4.05^ 

^•-'-'J-0.89V-2.00-' 
A qo+2.04/+3.72\ A 9q+1.55/-+2.96\ 
4-.3O_l.l6U2.i2J *-^'3_0.70Ul.50J 


WMAP7 + HST 
WMAP7 -^ HST + SDSS 
WMAP7 ■\- HST + CFHTLS 
WMAP7 + HST + SDSS + CFHTLS 
WMAP7 + HST + SDSS + SN + CFHTLS 


qq+l-76/+2.92x 
.3-yy_0.4ll-l.36J 
q Qq + 0.85/+1.78x 
0-yO_Q 73l,_i 42J 

q QO+1.04/+2.02-, , ,„+0.80/+1.62\ 
O-y0_g5i^_l 20J *-i'-0.67Ul.26J 
q Q-I+0.71/+1.38N q Q9+0.75/ + 1.42\ 
O-»J-_0. 68 V -1.20/ '3-^^-0.51 V- 1.03^ 
q Q9+0.61 /+1.33\ q Qi+0.67/ + 1.26\ 
"J-y^ -0.60 1-1.17 J '3-y J- -0.45 V -0.96 J 



C. Relativistic Species NcH 



In this subsection, we consider the constraints on the effective number of neutrino species, iVoff , assuming massless 
neutrinos. Since iVefi can be written in terms of flmh^ and the redshift of matter-radiation equality, Zs n, t here are 
strong degeneracies present between Ncs, the matter density, r2,„/i^ and the Hubble parameter Hq [3, l2J|. CMB 
constraints on Ncs can thus be strongly improved by combining them with measurements of the small-scale matter 
power spectrum such as those obtained from SDSS or CFHTLS. 

We find the WMAP7 data alone gives A^^eff > 2.84 at the 95% confidence level, consistent with the result derived 
by the WMAP7 team [01- Adding the CFHTLS data significantly improves the constraints on N^ff to 



A^cff = 4-38^ 



1.16 



/+3.72N 
1-2. I2J 



(18) 



at the 68% and 95% C.L. Figure [7] shows the two dimensional likelihood contours for both these constraints, with 
^max — 630, making it clear that the improvement is coming from degeneracy breaking- For comparison, when we use 



4-13 



+ 1.20^+4.05 
-0.89 



U2.00J 



the SDSS DR4 LRG power spectra data instead of CFHTLS, we obtain the very consistent A'cft 
(68 and 95% C-L-). 
WMAP7-hHST gives N^s = 3-99t.oAi{tlfe) (68 and 95% C-L.)- Adding the CFHTLS angular power spectra 



+0.85/ + 1.78\ 



tightens this to A^eff = 3-98Io;^i(lf;^o)- This is similar to WMAP7-HHST+SDSS, which gives N^s = '^■^^-orti 1A2) 



WMAP7-hHST-hSDSS-hCFHTLS gives A^cff = 3-91I^;^8(+i;^°)- Table |TTT] hsts these constraints on iVoff, while figure 
[5] shows one-dimensional marginalized likelihood constraints on A'off, i^mh^ and Hq, with ^max = 630- 
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FIG. 7: Marginalized one-dimensional and two-dimensional likelihood (1, 2 a contours) constraints on Qmh^ , Ho and A'^eff from 
WMAP7 (black solid lines) and WMAP7+CFHTLS (red dashed lines) data combinations for im^^ = 630. 
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FIG. 8: One-dimensional marginalized likelihood constraints on Ncs, ft-mh^ and Ho from different data combinations: 
WMAP7-(-HST (black solid hnes), WMAP7+HST-|^CFHTLS (blue dashed lines), WMAP7-HHST-HSDSS+CFHTLS (green 
dotted lines) and WMAP7-hHST-hSDSS-hSN-hCFHTLS (red dash-dot lines) for ^^ax = 630. 



WMAP7+HST+SDSS+CFHTLS+SN gives our most stringent constraint of 



iVcff = 3.92_„.6„ 



+0.61 (+1.33) (gg^^ ^^^ 95% c L^ 



(19) 



We also consider the ^max — 960 case and find that the constraint improves slightly: N^s = S-Q^1lo'%C^o'gQ) for the 
68% and 95% confidence levels. Our results are quite consistent with those of the ACT and SPT CMB experiments, 
which are A'cff = 5.3 ±1.3 and iVeff = 3.85 ±0.62 (68% confidence level), and, like them, display a slight preference for 
an extra relativistic relic [SlllS^. However, the standard value of N^s = 3.04 remains well within the 95% confidence 
intervals. 



VI. CONCLUSIONS AND DISCUSSIONS 



Measurements of the galaxy power spectrum can put strong constraints on the neutrino mass and effective number 
of neutrino species. Significant gains can be made by probing weakly non-linear scales and investigating redshift 
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evolution. We have used the improved parametrization developed by ref. yv\, who precisely calibrated the effect of 
massive neutrinos on the matter power spectrum using a suite of N-body simulations, to investigate non-linear scales. 

Although the non-linear matter power spectrum may be estimated accurately from simulations, connecting it to 
the galaxy clustering is challenging. In this study, we have assumed a simple bias model with no dependence on scale. 
To check the validity of this model, we fit the power spectrum with a simple two-parameter Q bias model in ACDM 
without massive neutrinos. We find that the best-fitting Q model has negligible scale-dependence to fc '--^ 0.6 h/Mpc. 
We stress that this is only a consistency check, but it demonstrates that a scale-dependent bias model is not needed 
to fit the data at the redshifts we consider, z > 0.5, even to weakly non-linear scales of fc ~ 0.6h/Mpc. These results 
are in part due to the range of galaxy luminosities represented in the CFHTLS sample. The CFHTLS data set thus 
probes an interesting regime for neutrino studies for z = 0.5 — 1.2 at small scales that can be modeled with confidence. 

We find that combining the CFHTLS data with WMAP already gives tight limits on the total neutrino mass. 
Adding additional low-redshift probes, including SDSS LRGs and SN, only marginally improves the constraints. 
We benefit considerably by extending the angular scale limit from ^max = 630 to 960. This extends the fit from 
k ^ 0.4h/Mpc to /c ~ 0.6h/Mpc. However, there is a strong correlation between galaxy bias and total neutrino mass, 
since both parameters modulate the amplitude of the galaxy power spectrum. Thus, future analyses can benefit from 
adding additional observables to constrain the galaxy bias. 

We further constrain the effective number of neutrino species. Again, CFHTLS complements the WMAP7 data by 
constraining fl„i leading to tighter limits on Ncs- We find that the constraining power of the CFHTLS data is similar 
to that of the SDSS LRG sample and the limits from both surveys are fully consistent. 

The main constraints derived in this paper are summarized as follows: CFHTLS, combined with WMAP7 and a 
prior on the Hubble constant provides an upper limit of Y^rriu < 0.29 eV and Ncs = 4.17^126 (^^ confidence levels). 
If we instead omit smaller scales which may be affected by non-linearities, these constraints relax to the following 
values: "^miy < 0.41 eV and A^eff = 3.98^1^20 C^^ confidence levels). By combining with the SDSS LRG matter 
power and SN luminosity distance moduli we further improve to ^ ttIj, < 0.27 eV and A'off = 3.91^ogg (2 a confidence 
levels) . 
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